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Abstract 


A method is presented for biasing spectral estimates to enhance detection of tonal signals 
against a background of broadband noise. In this method, a nonlinear average of an ensemble 
of individual spectral estimates is made where broadband noise energy is biased downward, 
pure tone energy is unbiased, and a mixture of the two is biased by an amount that depends 
on the ratio of tonal energy to broadband energy. The method is analyzed to provide 
estimates of the extent of tonal signal detection enhancement. 

Symbol List 


c 

contrast of adjacent spectral 

c q 

bands 

apparent contrast of adjacent 

Eq 

spectral bands exhibited by W q 
expectation operator for density 

Ei 

function of noise 
expectation operator for density 

function of noise plus tonal 

E2 

signal 

expectation operator for density 

/o 

function of tonal signal 
density function of single 

h 

spectral estimate of noise 
density function of single 


spectral estimate of noise plus 

h 

tonal signal 

density function of single 

Fo 

spectral estimate of tonal signal 
cumulative function of single 


spectral estimate of noise 



Fi 

cumulative function of single 
spectral estimate of noise plus 
tonal signal 

A 

density function of biased 
spectral estimate of noise 

ft 

density function of biased 
spectral estimate of noise plus 
tonal signal 

ft 

density function of biased 
spectral estimate of tonal signal 

A 

/o 

density function of ensemble 
averaged spectral estimate of 
noise 

/l 

density function of ensemble 
averaged spectral estimate of 
noise plus tonal signal 

ft 

density function of ensemble 
averaged spectral estimate of 
tonal signal 

1*1 

degenerate hypcrgeometric 
function 

G 

arbitrary function 

I v 

modified Bessel function 

m 

number of spectral estimates in 
an ensemble 

P(D) 

probability of detection 

P(FA) 

probability of false alarm 

9 

power of a power function 
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R 

Rq 

R 

R 

T 

T 

T 

W 




W\ 


signal-to-noise ratio 

apparent signal-to-noise ratio 
exhibited by W q 
apparent detection signal-to- 
noise ratio from biased spectral 
estimate processing gain 
apparent detection signal- 
to-noise ratio from ensemble 
average processing gain 
detection threshold 

detection threshold of biased 
spectral estimate 
detection threshold of ensemble 
averaged spectral estimate 
statistic composed of an inverse 
function of the average of 
an arbitrary function of an 
ensemble of single spectral 
estimates 

statistic composed of the q ^ 1 
root of the average of the power 
function, of power q , of an 
ensemble of single spectral 
estimates 

statistic composed of the 
arithmetic mean of an ensemble 
of single spectral estimates 
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Wq 

statistic composed of the 
geometric mean of an ensemble 


of single spectral estimates 

W - 1 

statistic composed of the 
harmonic mean of an ensemble 
of single spectral estimates 

W-oo 

the limit of W q as q — > — oo, 
equivalent to the first order 
statistic of an ensemble of single 
spectral estimates 

X 

spectral estimate consisting of 
the power spectrum, or squared 
magnitude of the Fourier 
transform, at a single frequency 

x i 

t 1 ^ spectral estimate of an 
ensemble 

*(1) 

first order statistic 

*m i 

ensemble of m independent 
single spectral estimates in a 
frequency band 

x (m) 

ensemble of m ordered sin- 
gle spectral estimates in a fre- 
quency band 

7 

Euler’s constant 

r 

gamma function 

6 

Dirac delta function 

C 

dummy variable of integration 

V 

dummy variable of integration 
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root of false alarm equation 
noise power 
tonal signal power 

Introduction 

Spectral estimates for helicopter acoustic signatures are generally made by calculating a 
set (or ensemble) of presumably independent power spectra using a finite Fourier transform 
and then linearly averaging those spectra. By averaging, the uncertainty of the spectral 
estimate is reduced, and a peak (tone) in the spectrum is made more distinguishable from 
the random background noise. Although the process of averaging individual spectra is linear 
(in terms of squared pressure), the individual spectra and the average spectrum are biased 
in the region of a tone. Despite the bias associated with a finite Fourier transform estimate 
of a spectrum, the process of linearly averaging a set of spectra introduces no further bias. 
Also, by reducing uncertainty, the average provides a better estimate of spectral levels than 
an individual spectrum. 

By assuming that a random process is essentially broadband in nature (and avoiding 
the bias associated with a finite Fourier transform by letting the integration time approach 
infinity), it is possible to derive the probability density function of a single spectral estimate 

/„(,)= 

where x is the spectral estimate (in terms of squared pressure) in a frequency band, and 
a\ is the portion of broadband (noise) power in that band. It is also possible to derive the 
probability density function of a single spectral estimate of a composite process composed 
of both sinusoidal (tonal) and broadband (noise) terms 
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where <r a = P 2 / 2 is the power of a sinusoidal signal of pressure amplitude P, and Iq() is the 
zero order modified Bessel function [1], When the signal-to-noise ratio becomes very large, 
the probability density function approaches that of a sinusoidal signal without noise 


f 2 (x) = 6 (x- afj 


where £() is the Dirac delta function. 

The different functional forms of the noise probability density and the signal-plus- 
noise probability density suggest that it may be possible to enhance the difference in their 
statistical behavior to provide a better apparent signal-to-noise ratio at the expense of 
generating biased spectral estimates. Specifically, an estimation method that tends to 
emphasize smaller values of an ensemble of individual spectra, rather than equally weighting 
all spectra in a spectral average, should bias a pure noise spectrum more than one composed 
of both tonal signal and noise. Ideally, the estimate should be completely unbiased for a 
sinusoidal signal with no noise. 

A Nonlinear Estimate Method 


A direct method to form a biased spectral estimate is to evaluate the inverse function of 
an averaged function of the individual spectra 


W = G~ l 



where X{ is the member of an ensemble of m independent spectral estimates in a particular 
frequency band 


G 



Xm — {*1) x 2> • • • > 


A power function 


G(x) — x q 


where q < 1, forms a relatively simple statistic 


IV, = W,(Xm) = (iy'x? 


m 


m \ , 
*=1 


with a very useful characteristic, namely, that W qi < Wq 2 if 91 < 92 an( i the X{ are not 
all equal. For q = 1, W\ is the arithmetic (unbiased) mean, as q -* 0, Wq approaches 
the geometric mean, and for q = -1, W-\ is the harmonic mean. Each frequency band 
in an ensemble of spectra is processed independently of the others so that the collection of 
frequency bands so averaged forms an ensemble averaged spectral estimate. 


Analysis of the Method 

Unbiased Estimate 

The expected value of a linear ensemble average (unbiased) spectral estimate for broad- 
band noise alone is 

f°° 2 

E 0 [Wi] = E 0 [x}= xf 0 (x)dx = *l 
\ JO 

For signal plus noise, the expected value is 
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roo 

E\[W\) = Ei[x)= / xfi(x)dx = erj + trj 
JO 

and for signal alone, the expected value is 

f°° 

^2(W^l] = ^2l*l = / xf2(x)dx = a* 

JO 

Indeed; for signal alone, the joint probability density function of m independent spectra is 
given by 

m 

/2(%) = n 8 (*» ~ & >) 

1=1 

so that the expected value for signal alone of any estimate Wq is given by 

W q (x m ) / 2 (x m ) dx m = 

which indicates that W q is, as desired, an unbiased estimator for any value of q when a signal 
is present without noise. 

Geometric Mean 

To determine the effect of the biased estimator, Wq , on the signal-to-noise ratio when 
q ^ 1, the expected values of the estimator for noise and signal-plus-noise must be compared 
with the expected values of an unbiased estimator, W\. The simplest analytical case is q = 0. 
Because Wq is the geometric mean, it can be expressed as 

Wo(x m J = •••Zm 

so that the expected value for noise only is 
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BoIWo] = {Bo [ VS]} m = <^[r(l + 


m 


and the expected value for signal plus noise is 


Bi[Wb) = {£i[VS]} m 

where the signal-to-noise ratio is R = tr^/o^, and the degenerate hypergeometric function 


16 


1^1 (“, P, z ) = 1 + 




which must be evaluated numerically. 

The signal-rto-noise ratio of an unbiased estimate can be expressed as 


_ filial , d 

BolH'l] cl 


so the apparent signal-to-noise ratio of a biased estimate should be 


, B|[W„1 

1 BolW,] 


For 9 = 0, the apparent signal-to-noise ratio is then 



Ei[W 0 ) 

Eq[W 0 } 




-1 


which depends only on the unbiased signal— to— noise ratio, R t and the number of independent 
spectra, m, included in the estimate. 
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Harmonic Mean 


A more involved analysis is required for the case of q = —1. For noise alone, the joint 
probability density function of m independent spectra is given by 




The expected value for noise alone of an estimate is given by 


E 0 [W-i] = f W_i(x m )/o(x m )dx m 
JO 

Substituting the appropriate expressions in the equation above gives 


Eo[^-i] = 


mal 


r nJ Q z i 1 j e dz\dz2...dzm 

Only for the simplest non-trivial case, m = 2, can a closed form solution be derived 

E 0 [H'.,] = 2 Q) 

No attempt was made to analyze signal-plus-noise. 


First Order Statistic 


If an ensemble of independent spectral estimates in a specified frequency band is given by 


x m — {*1 j x 2j • • ■ t ^m} 


then, when these spectral estimates are arranged in ascending order, the (dependent) ordered 
samples are given by 
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where 



*(1) ^ *(2) ^ ^ x (m) 

and x^j is referred to as the order statistic. It can be shown that 

W-oo = lim W q = z(i) j 

which means that the smallest possible value of the biased estimate is the first order statistic 
or the smallest of the spectral levels in a frequency band. 

Because the first order statistic, *(1)1 provides a lower bound for the suggested biased 
estimator, an analysis of the apparent signal— to-noise ratio for this estimator will provide 
an indication of the maximum signal-to-noise ratio enhancement. The probability density 
function of the first order statistic of m independent samples is given by 

/(*(!)) = m[i-F(x = X(i))] /(x 

where F(x) is the cumulative distribution function 

F(x) = f /«K 

Jo 

of a single spectral estimate [ 2 ]. For noise alone this gives 

Fq{x) = 1 - e”*/*" 

which, in turn, gives 
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which has an expected value of 


roo 

•Eo [*(i)] = J Q v'MvWv = *1/™ 


For signal-plus-noise, the cumulative distribution function is 


0i(*)= f /i(CK 

Jo 


which, in turn, gives 


/l(*(i)) = m 1 “ J Q ( } MO d C fl( x = *(!)) 


which can be integrated numerically to determine the expected value 


0l[*(l)] = m J q 


-i m—1 


fi(v) d v 


The apparent signal-to-noise ratio of the biased estimate provided by the first order 
statistic can then be expressed as 


n SllH'-Oo] ■ , 

“ £o|IV-«] Bo[*(l)] 

_2 roo r pr\ I m—1 

= — / vU~ fi(C) d C f\{n) d i - 1 
°n JO L JO J 

which gives the best signal-to-noise ratio enhancement achievable by these methods. When 
the unbiased signal-to-noise ratio is very large, the apparent signal-to-noise ratio is bounded 





E 2 [W- oq] 
Eo[W-oo] 



— 1 = mR — 1 


On a logarithmic scale, the best signal-to-noise ratio improvement that can be expected is 
given by 


AS NR = 10log 10 (i2-oo/#) « 10log 10 (m)dB 


Not surprisingly, the method of biased spectral estimates for enhancing signal-to-noise ratio 

works best when the unbiased signal-to-noise ratio is very high. 


Signal— to— Noise Ratio and Contrast 

A useful concept for examining the effect of signal-to-noise ratio enhancement is contrast. 
If two adjacent frequency bands contain the same level of broadband noise while only one 
contains a tonal signal, then the signal-to-noise ratio in classical analysis is given by 


_ Ell*] 

E 0 [x] 




The contrast can be expressed as the ratio of the two spectral levels (or the difference between 
the two spectral levels on a dB scale) 


c=|^= i +4= i+fi 

E 0 [x] °n 

For the case of a biased estimate, the same basic definitions hold true. The signal-to-noise 
ratio is given by 


_ A [*(!)] x 
£ol*(i)] 
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while the contrast is given by the ratio of the two spectral levels 


-OO 


= £ihyi i +iJ _ 


For example, when the unbiased signal-to-noise ratio is 0 dB, R - 1, then the unbiased 


contrast is about 3 dB, (7 — 2. 

The utility of biased spectral estimates can be shown by comparing the apparent signal- 
to-noise ratio of a biased estimate with that of a classical one, R q /R, or by comparing the 
apparent contrast of a biased estimate with that of a classical one, C q /C. The only two 
parameters analyzed were q = 0 and q = -oo. Figure 1 shows the signal-to-noise ratio 
enhancement in dB that can be expected from using the harmonic mean, W 0 , instead of 
the arithmetic mean, W X , for several different signal-to-noise ratios. It is clear that the 
greater the signal-to-noise ratio, the greater the enhancement. However, regardless of the 
signal-to-noise ratio, the limit to the extent of enhancement is given by 


lim 

R—*oo 


Rp 

R 



- It is also clear that the greater the number of spectra included in the technique, the greater 
the enhancement. However, regardless of the number of spectra, the limit to the extent of 

enhancement is given by 


lim % = el a 2.51dB 
m— *oo R 

where 7 is the Euler constant. Figure 2 shows the corresponding contrast enhancement. 

Figure 3 shows the signal-to-noise ratio enhancement in dB that can be expected from 
using the first order statistic, * (1) = W-oo, instead of the arithmetic mean, W Xt for 
several different signal-to-noise ratios. The greater the signal-to-noise ratio, the greater 
the enhancement, but the limit to the extent of enhancement regardless of the signal-to- 


noise ratio 



R — oo 

lim — — = m 
R—*oo it 

Figure 4 shows the corresponding contrast enhancement. 

Signal Detection and the Threshold Effect 

A simple analysis of the first order statistic method shows that there is a signal— to - noise 
ratio enhancement based on the mean value of the distribution of the statistic. However, 
signal detection depends on exceeding a threshold which is determined from an acceptable 
probability of false alarm, P(FA). Because the detection threshold depends both on the 
mean and the variance of the statistic, calculating a signal— to— noise ratio enhancement for 
detection analysis that depends only on the mean can be misleading. For a single spectral 
estimate, the noise density function is given by 



from which the threshold, T, can be determined by solving 


so that 


P{FA) = r fo{x)d 
JT 



r = -4in|P(^)l 

The probability of detection when a sinusoidal signal is present can then be written 

l 

fl(x)dx 

where the density function of signal plus noise is 


P{D) = 
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/I(I) = ( 4 ) 10 (^f ) 

The noise density function for the first order statistic, from an ensemble of m spectra 

is given by 

/«(*<!)) = Q) 

The threshold, T, can be determined by solving 

- 

p ( FA )= . M x (l)) dx ( 1) 

•/T 

which gives 


t = -<rl\n[P(FA)]/m 


The probability of detection when a sinusoidal signal is present can then be written 


r °° . 

p ( D ) = /. /i( I (i)) da; (i) 

J T 


where the density function of signal plus noise is 


/l(*(l)) = m 1 “ J Q 1 } fl( x = *(!)) 


For the statistic W\ t which is a linear average of an ensemble of m spectra, the noise density 
function is given by 


/o(Wi) = 


m fmWi\ m 1 e -mW 1 /al 

ff n r ( m )Un / 
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Similarly, the threshold, T, can be determined by solving 


P(FA) = /°° fo(Wi)dWi 
JT 


so 


that 


T = (</m 


where £ is found by solving 


m-i £jb 


P(FA) = e-« £ 


k=0 


The probability of detection when a sinusoidal signal is present can then be written 


P(D) = /°° fi(W\)dWi 
JT 

where the density function of signal plus noise is 

= 

The simplest way to examine the effect on detection of using or W— oo, rather than 

W\ t is to assume that the estimates are approximately unbiased for a sinusoidal signal plus 
noise, as when the signal-to-noise ratio is very high, and examine the ratio of the respective 

thresholds 



f/T=-l n [P{FA)}/i 

where £ is as above. Figure 5 shows the ratio of the classical threshold to the first order 
statistic threshold, in dB, for various probabilities of false alarm. Because these curves are 
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only appropriate for high sign— to— noise ratios and detection is desired for marginal signal 
to-noise ratios, this approach can be misleading. 

A less extreme simplification is to assume that the W—oo estimate is biased when a 
sinusoidal signal is present and that detection occurs when the expected value of signal plus 

i , 

noise is equal to the threshold. For W\, detection would then occur when 

T = a 2 +cr 2 

from which the classical signal-to-noise ratio can be derived 

. R — i/m — 1 

where i is, again, as above. For W— oo» detection would occur when 


T — 





m— 1 

fl{x)dx 


and the amount the biased technique would enhance detection is then given by the ratio 
of the two signal— to— noise ratios. Determination of the signal— to— noise ratio, R, requires 
finding the root of the equation 


where 


0 = In [P(FA)\ 

+ m 2 a{R,() d ( 


m— 1 

g(R,v) d V 


g(R, V ) = e-b+V 

The signal-to-noise ratio enhancement is then R/R. Figure 6 shows curves of the signal- 
to-noise enhancement for detection purposes. Because the first order statistic shows greater 
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variability than the classical average despite the fact that is has a much lower mean, the 
enhancement shown is not as great as indicated by simple signal-to-noise ratio or contrast 
methods. 

Although a probability of detection was not specified, using the expected value of the 
density function to define detection is not entirely unfounded. For a symmetric density 
function, the expected value would yield a probability of detection of 50%. The classical, 
W\, density function is asymmetric but can be shown to approach a symmetric Gaussian 
form when either the signal-to-noise ratio is high or when the number of spectra in the 
average is great. There should, however, be minor differences between the detection signal- 
to-noise ratio enhancement curves shown here and those curves which might be derived for 
a fixed detection probability. 

Conclusions 

A method was presented for calculating biased spectral estimates that enhance tonal 
signals against a background of broadband noise. The method was shown to differentially 
bias different mixtures of broadband noise and pure tones as a function of the ratio of tonal 
energy to broadband energy. The method was analyzed and shown to provide the best 
enhancements for large ratios of tonal energy to broadband energy. The method provides 
some enhancement for any tone, but is unable to significantly improve detection of tones 
that are very much lower than broadband noise levels. 
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S/N Ratio Enhancement, dB 



Figure 1. Signal-to-noise ratio enhancement from the harmonic mean. 


20 



Contrast enhancement, dB 



Figure 2. Contrast enhancement from the harmonic mean. 
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S/N ratio enhancement, dB 



Number of samples 

Figure 3. Signal-to-noise ratio enhancement from the first order statistic. 
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Figure 4. Contrast enhancement from the first order statistic. 
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Figure 5. Detection threshold enhancement from the first order statistic. 
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Figure 6. Detection signal-to-noise ratio enhancement from the first order statistic. 
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